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Abstract 

We investigate the case of a homogeneous tachyon field coupled to gravity in a spatially 
flat Friedman-Robertson- Walker spacetime. Assuming the field evolution to be exponentially 
decaying with time we solve the field equations and show that, under certain conditions, the 
scale factor represents an accelerating universe, following a phase of decelerated expansion. We 
make use of a model of dark energy (with p = —p) and dark matter (p = 0) where a single scalar 
field (tachyon) governs the dynamics of both the dark components. We show that this model 
fits the current supernova data as well as the canonical ACDM model. We give the bounds on 
the parameters allowed by the current data. 
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I. INTRODUCTION 

n 

The observational fact that the present universe is accelerating |l| has led to several 
theoretical speculations about the nature of the acceleration and the cause behind it. 
No completely consistent model exists till date, though a plethora of them with intrinsic 
advantages and discrepancies have been proposed (for a sample of these in recent literature 
see Q). The cause of the acceleration is attributed to the presence of a mysterious entity 
named dark energy which dominates (is about 70 percent) the matter-energy content of 
the present day universe jjj. In order to drive the acceleration the Dark energy needs 
negative pressure. For an equation of state of the form p = wp observations place bounds 
on w : —1.33 < w < —0.8. Thus models respecting p = wp must be consistent with 
these bounds. Among many, the simplest is Einstein's old idea of an universe with a 
cosmological constant, for which w = — 1. In a more general scenario the case of a time 
dependent w(t) has also been looked at to some extent. 

The majority of dark energy models involve the temporal dynamics of one or more 
scalar fields with varied actions and with or without a scalar potential 4]. A sub-class 
of these scalar field models includes exotic models where the kinetic energy term may 
become negative (the phantom) .5] or the scalar action may be of a non-standard form 
(the tachyon condensate). The latter arises in the context of theories of unification such 
as superstring theory p. In this article, we shall concentrate on the the cosmological 
consequences of the class of models where the scalar tachyon condensate appears in the 
matter part of the Einstein equations. 

Let us first briefly review the origins of the effective action for the tachyon condensate. 
It is known from the early days of string theory that the spectrum of strings do contain 
a negative mass squared vibration called the tachyon. For strings attached to Dirichlet 
(D) branes such tachyonic modes reflect D-brane instability. The low energy effective 
action for the scalar tachyon field around the tachyon vacuum thus contains a tachyon 
potential which has an unstable maximum resulting in a negative mass square value in the 
mass term. The scalar tachyon rolls down from the potential maximum to a more stable 
configuration thereby acquiring a stable vacuum expectation value. This is the tachyon 
condensate. The questions in cosmology are : (i) Can such a rolling tachyon drive the 
acceleration of the universe? (ii) Can the tachyon condensate be of use in describing all 
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that is described in cosmology using scalar fields (inflation, dark matter etc)? A distinct 
advantage of using the scalar tachyon to understand cosmology is the fact that we do not 
need to answer the embarrassing question 'where does the scalar field come from?' Note 
that this is reminiscent of the oft-quoted advantage of string (dilaton) cosmology too Q| • 

There have been several articles in the literature where the cosmological consequences 

an 
. It was shown in [9] that the effective 

tachyon action can lead to an energy-momentum equivalent to a gas with zero pressure 

(dust) but nonzero ener gy density (modulo some assumptions and re-definitions). It 

has also been noticed in that, modulo some assumptions, the tachyon field can also 

be made use of in describing dark matter at galactic scales. In addition, the energy 

momentum tensor of the tachyon condensate can be split into two parts-one with zero 

pressure (dark matter) and another with p = —p. This facilitates the description of dark 

energy and dark matter using a single scalar field-a fact we shall test in detail for our 

model in this article. Tachyons and their dynamics in the context of string theory have 

been reviewed in 

Most of the work on tachyon cosmology assumes a tachyon field varying linearly with 
time. In addition, the choice of the tachyon potential is also assumed to be either expo- 
nential (decaying) or inverse quadratic Alternative potentials satisfying the general 
features (derived from string field theory) have also been proposed and analyzed [121 ] . In 
our work, we shall assume the tachyon field to be decaying exponentially with time. We 
will see that the scale factor can have an accelerating behavior for certain values of the 
constants which appear in our solution. The tachyon potential for such scenario does 
satisfy the general requirements. In addition, we find different types of solutions with 
the exponentially decaying tachyon field. Finally, we use Supernova data to constrain the 
parameters that appear in our solution. 

We choose to work with units c = 1 though, wherever necessary we go back to the 
actual units and dimensions. 



II. THE ANALYSIS 



The action for the homogeneous tachyon condensate coupled to gravity is known to be 
of the form j^| 
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where k = 8irG and V(T) is the tachyon potential. We choose the tachyon field to be 
time dependent T = T(t) in a spatially flat 3+1 dimensional FRW background. The line 
element is 

ds 2 = -dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) , (2) 

where a(t) is the scale factor. The exact form of the potential is not known. However as we 
argued above, (i) it should have an unstable maximum at the origin and (ii) should decay 
exponentially to zero as the field goes to infinity. Moreover, the slope of the potential 
should be negative for T > Q|. Assuming the field stress energy tensor of the form, 
diag (p,p,p,p), the energy density and pressure obtained from eqn (JTJ) are as follows: 

V(T) 
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Ignoring the cosmological constant the two independent components of the Einstein 
field equation give: 



kV(T) 
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For the three unknown functions (V(T), a(t) and T(t)) we have two equations and to solve 
for these we need to choose one of them suitably. Since no initial or boundary conditions 
are imposed, some arbitrary parameters are expected to appear in the solutions but we 
will see that these will not be completely arbitrary in order to give consistent, physically 
acceptable solutions. Also these parameters will be finally estimated by making a fit with 
observed supernova data. Depending on one's choice of field or potential or the scale 
factor, a large class of consistent solutions can be constructed, (e.g. 3 discusses exact 
solutions starting with chosen form of inverse square potential, and [l^J discusses solutions 



based on choice of T(t), whereas |15( discusses solutions obtained from a different choice 
of the potential (exponentially decaying)). One problem that might arise in the approach 
we adopt, is that the resulting potential may not be physical. This will be revealed as we 
proceed with the calculations. Assuming the tachyon field to be decaying exponentially 
with time we have 

T(t) = T e~ at , (7) 

where a is a positive constant and T Q is the value of the field at t = 0. Using this form of 
the field in ©,© and eliminating V(T) one gets the Hubble parameter for the universe: 

i = 1 (R ) 

a C _ ^Il e -2at ' { } 

where C is an arbitrary integration constant. Incidentally, a physically acceptable theory 
is obtained only when C is positive. (C < gives collapsing universes.) Integrating the 
last relation and using the fact that a(t = t , present epoch) = 1 we get 
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where t is the age of the universe in the model under discussion. To make calculations 
simpler we now define three dimensionless parameters: 

> = (10) 

Using these one can write 

a(t)=( e —P) S . (11) 
\ r — p J 

The other unknown function, the tachyon potential can be obtained from 



V (T) = t c _ ¥r / )3 • (12) 



The potential conforms to the nature of other tachyon potentials provided the maximum 
of potential occurs at a higher value of the field. Note that in our case the tachyon field 
decays with time. This requirement leads to the condition Ca < 3. In terms of the 
dimensionless parameters the time dependence of potential can be expressed as 

3(1 _ 2p e -2ani/2 

v{t) = i {1 -V J)» • (13) 



The potential will be a real function for all values of time if || < 1 i.e. p < y. 
However, from the expression for the scale factor we note that a(t = 0) = (big bang) 
would imply a choice p = 1 (and hence s > 2/3). These two conditions give the required 
potential for our purpose (see fig ll|2j) which like the 1/T 2 potential is singular at the 
origin (in time). It can be easily seen that the isotropic pressure and energy density also 
remain well behaved under these conditions. Another important feature of the theory is 
the deceleration parameter which is given by 
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and represents an accelerating universe at late times 
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Field T in units of \ / a* 



FIG. 1: We plot the potential V(T) as a function of the field T in arbitrary units. The values 
of the model parameters for this and the next figure are taken from best fit to the current 
supernova data given later. 



III. FITTING THE SUPERNOVA DATA : PARAMETER ESTIMATION 

The tachyon condensate has been proposed as a candidate for dark energy. The above 
exact solution, however, does not assume anything other than the tachyon condensate as 
a source in the Einstein equations. Therefore, if the above scale factor is assumed to be 
correct then we have to believe in a universe where dark energy is the sole matter source 
(which is unrealistic). To remedy this deficiency we can propose the following. 
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FIG. 2: Potential as a function of time. The field rolls down from top of the potential and after 
infinite time interval reaches a constant value. 

(a) Add other matter fields (dark matter + A) and obtain the scale factor analyti- 
cally/numerically. Then obtain the unknown parameters by fitting with supernova data. 
Finally, use these parameters to ascertain the proportions of dark matter and dark energy 

(b) Use the split of the energy density and pressures of the tachyon condensate into 
dark matter (pressure- less) and dark energy (with w — — 1) parts. This goes as follows : 



P — Pdm + Pde 



P = Pdm + Pde , 



(15) 



where 
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In this way one should be able include the contributions of both dark matter and dark 
energy. 

We shall first explore the second possibility. Using the relation between redshift and 
scale factor, 1 + z = a /a(t) = l/a(t) one gets time-redshift relation in the present case 
as: 



t(z) = In 
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We notice one important restriction following from the last relation. If we impose the 
condition that as t — > 0, z — > oo then we must make p = 1. We will use this value when 
calculating some parameters of the model of the universe later. From this, one can express 
Qdm and VLde as function of redshift as follows: 



^dm(z) = f 2 = Tlo? exp (-2at(z)) (19) 
n DE {z) = l-n DM {z) (20) 

where t(z) is as given in EqUHl Notice that in (fTH|) and (J2U|) to denotes (present) age of 
the universe. We plot Qdm{z) and Q,de{z) in fig (jHJ). 
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FIG. 3: This plot shows how Qdm and £Ide evolve with redshift in our best fit model. The 
universe is completely dominated by matter beyond a redshift z ~ 10. 

One way to test the validity of our model is to calculate its prediction of &dm and 
compare with the known value (in the LCDM model). It should be noted that in this 
model these energy densities are only 'effective' energy densities, and these components 
do not truly behave like the usual dark matter and dark energy. However, since we can 
think of these densities as the appropriate terms in the Taylor expansion of the Hubble 
parameter in our model, it is likely that the local behaviour that fits the supernova data 
would be similar, and hence these densities might compare well. Another constraint that 
our model must satisfy is the age of the universe, which is known from independent 
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measurements. One can define the so-called Hubble free luminosity di distance using 
Dl = cH^di- Using the definition of luminosity distance 

/ H 



d L z) 



:i + z) f 

Jo 



dz' 



H(z' 



21) 



we derive the expression for the same in our model as: 
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In order to test the suitability of the present model of a tachyon condensate as the origin 
of dark matter and dark energy we have used the observational Gold dataset comprising 
157 Type la supernovae The Hubble constant free distance modulus is given by 



//(*) = 5 log 10 (rf L (z)/Mpc) + 25 . 

For model fitting we need to consider the likelihood function given by 

X 2 (H ,p/r,s) 
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where /ij and Z{ are the observed distance modulus and the redshift, and /ifi t is the ex- 
pression Eq |2H1 with the luminosity distance given in Eq 21. The posterior probability for 
the model parameters (H ,p/r, s) is given by 



P(iWr,s)ocexp(-|-)Pr(#o), 



(26) 



where Pr(ifo) refers to the prior on the Hubble parameter. For our analysis we have used 



the Gaussian prior, H = 66 ± 6 (|l2|)- We need to define the bounds for the three- 
dimensional volume in (H ,p/r, s). The lower and upper bounds on H are taken at 50 
and 80 . It is obvious from the following expression that p/r < 1. 

1 



(27) 



C(l-p/r) ' 

But for q there is no upper bound. So we choose a sufficient range for these parameters 
(0.1 < p/r < 0.91 and 0.66 < s) in such a way that p/r x s is of the order unity. 

Starting from the best-fit p/r, s (the value at the maximum of the probability distribu- 
tion), we may move down till 68% probability is enclosed under the surface and obtain the 



la bound on the parameters. The 2a and 3a bounds can be similarly obtained. In Fig0] 
we plot the 68 % (solid), 90 % (broken) and 95 % (dotted) contours. The maximum of the 
likelihood surface is at s = 0.66 and p/r = 0.385. The Xmm/d-o.f. for this fit is 1.16. We 
simultaneously obtain Hq = 64.00. One can see that the contours are not closed in the 
allowed region in the parameter space. For this fit, The value of Ca = 0.76 which satisfies 
the constraint stated earlier. In Fig El we compare the present fit with the canonical best 
fit AC DM model. 
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FIG. 4: The innermost, middle and the outermost contours correspond to 68%, 90% and 95% 
confidence level respectively. Peak of likelihood surface corresponds to p/r=0.385, s=0.66. The 
asterisk marks the minimum of x 2 - As discussed in the text, to keep the potential physical at 
all times we must impose the constraint s > 2/3. We find that the maximum of the probability 
density lies at the edge of this interval. 



Accepting the values of the parameters obtained from fit we can now estimate various 
properties of the universe and compare with known values. Since best fit p/r = 0.385, s = 
0.66, H = 64.0 and we argued p — 1, it follows, r = 2.6 and a = 0A6H . Thus 
the age of the universe according to this model is t = log(r)/2a = 15.2 Gyr. Also 
(floM)o — 0.39, (floE)o — 0.61, q = p/rs — 1 = —0.42. These values are close to the 
currently accepted values. The deceleration parameter and redshift in this model are 
related by 

0.58 , . 

q(z) = ; t^t 1 , 28 

qK ' 0.615(1 + z)~ 1 - 5 + 0.385 ' V ; 
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FIG. 5: The green curve is the Hubble plot for the best fit model with p/r = 0.385, s = 0.66 
and Ho = 64.0. For comparison we also show the best fit AC DM model (blue curve). 

and is found to represent an universe which is accelerating at present epoch but was 
decelerating in remote past (see fig (JHJ)). This fact is also in agreement with current 
picture. The transition takes place at z — 1.15. Note also that apart from the best-fit 
values there are a whole range of values of the parameters which are allowed at different 
confidence levels. 
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FIG. 6: The deceleration parameter calculated for the best fit model as in previous figure. We 
see that the universe starts decelerating around z ~ 1. 

It is obvious that there are many inadequacies in this model. For instance in the remote 
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past there should be a radiation-dominated phase and the amount of radiation should 
decay to very small values in the present epoch. At best, this model can accomodate 
the late-time accelerated phase which, by and large, has been our goal. We have an 
analytic expression for the scale factor and the tachyon field. In addition, we have made 
a simplistic assumption about dark matter and dark energy being born out of the same 
scalar field. In some sense, our model tests how good this assumption is. The fact that 
the model fits the data quite well is appealing but to improve on it we must add further 
necessary features which are required for a realistic cosmological model. Finally, it might 
be worthwhile looking at cosmological perturbations in this model and try to use WMAP 
data to set more stringent bounds on the parameters [tt! 

Acknowledgments 

Work of AD and SG is supported by a Junior Research Fellowship by the Council of 
Scientific and Industrial Research, India. 



[1] A. G. Riess et al, Astron. J. 116, 1009 (1998) [ArXiv: astro-ph/9805201 ; A. V. Filippenko, 
A. G. Riess, Phys. Rept. 307, 31-44 (1998) [ArXiv: astro-ph/9807008 

[2] V. Balasubramanian, Class. Quant. Grav. 21, S1337 (2004) [ArXiv: |hep-th/0404075] ; 
E. S. Santini, G. A. Lemarchand [ArXiv: |astro- ph/0410056|; R. Wi gmans, [ArXiv: 
|astro-ph/040"9033"l ; A. B. Balakin, Gen.Rel.Grav. 36, 1513 (2004) [ArXiv: |gr-qc/0409024] ; 



S. Ray, U. Mukhopadhyay [ArXiv: astro-ph/0407295 
[3] C. Armendariz-Picon, V. Mukhanov, P. J. Steinhardt, Phys.Rev.Lett. 85, 4438 (2000), 
[ArXiv: |astro-ph/0004134| ]; V. Sahni, A. Starobinsky, Int.J.Mod.Phys. D9 373-444 (2000); 



P. J. E. Peebles, B. Ratra, Rev.Mod.Phys. 75, 559 (2003) [Arxiv: astro-ph/0207347 ; T. 

Padmanabhan, Curr. Sci.88, 1057 (2005) 
[4] I. Maor, R. Brustein, Phys.Rev. D67, 103508 (2003) [ArXiv: hep-ph/ 0209203 ; V. H. 

Cardenas, S. D. Campo, Phys.Rev. D69, 083508 (2004), [ArXiv: astro-ph/040"l03"T] ; P. G. 

Ferreira, M. Joyce, Phys. Rev. D58, 023503 
[5] R. R. Caldwell, M. Kamionkowski, N. N. Weinberg, Phys.Rev.Lett. 91, 071301 (2003) 



[ArXiv: |astro-ph/0302506] ; V. B. Johri, Phys.Rev. D70 041303 (2004), [ArXiv: 
|astro-ph/031 1293| ; V. Sahni, [ArXiv: |astro-ph/0403324| ; S. Nojiri and S. D. Odintsov, 
Phys. Letts. B562,147 (2003); S. Nojiri and S. D. Odintsov, Phys. Lett. B571 (2003) 1 



[6] A. Sen, ArXiv: |hep-th /0410103 



[7] A. A. Tseytlin, String cosmology and dilaton, in String quantum gravity and physics 
at the Planck scale, Erice, 21-28 June 1992, ed. N. Sanchez, World Scientific, 1993, 
[ArXiv: hep-th/9206067 ; J. E. Lidsey, D. Wands, E. J. Copeland, Superstring Cosmol- 
ogy, Phys.Rept. 337, 343-492 (2000) [ArXiv: |hep-th/ 9909061] 



[8] G. W. Gibbons, Phys.Lett. B537 1 (2002) |ArXiv:hep -th/0204008 ; G. W. Gib- 
bons, Class.Quant.Grav. 20, S321 (2003) [ArXiv: |hep-th/030lTT7| ; P. F. Gonzalez- 
Diaz, [ArXiv: |hep-th/0408225] ; B. C. Paul, M. Sami, Phys.Rev. D70, 027301 
(2004) [ArXiv: |hep-th/ 03 12081 M. Fairbairn, M. H.G. Tytgat, Phys.Lett. B546, 1 
(2002) [ArXiv: |hep-th/0 204070 ; G. Calcagni, Phys.Rev. D69, 103508 (2004), [ArXiv: 
hep-ph/0402126 ; G. Felder, L. Kofman, A. Starobinsky, JHEP 0209, 026 (2002) [ArXiv: 



hep-th/0208019 ; A. Mazumdar, S. Panda, A. Prez-Lorenzana, Nucl.Phys. B614, 101 
(2001) |ArXiv:hep-ph/0107058| ; X. Li, D. Liu, J. Hao, [ArXiv: |hep-th/0207146l 
[9] A. Sen, JHEP 0207 065 (2002) [ArXiv:hep-th/0203265] 
[10] T. Padmanabhan and T. Roy Choudhury, Phys.Rev. D66 081301 (2002) [ArXiv: 
|hep-th/0205055| ] 

[11] T.Padmanabhan, Phys.Rev. D66, 021301 (2002), [ArXiv: |hep-th/0204"T50| ; J. S. Bagla, H. 

K. Jassal, T. Padmanabhan, Phys.Rev. D67, 063504 (2003) [ArXiv: |astro-ph/0212198| , 
[12] C. Kim, H. B. Kim, Y. Kim, Phys.Lett. B552 111 (2003) [ArXiv: |hep-th/021oToT| ; 
[13] D. A. Steer, F. Verizzi, Phys.Rev. D70 043527 (2004) [ArXiv: |hep-th/0310"l39"l 
[14] A. Feinstein, Phys. Rev. D66, 063511 (2002) 

[15] M. Sami, P. Chingangbam, T. Qureshi, Phys.Rev. D66 043530 (2002) [ArXiv: 
|hep-th/0205179 



[16] A. G. Riess et. al, Astrophys. J. 607, 665 (2004) 
[17] M. Tegmark, et al., Phys.Rev. D69, 103501 (2004) 

[18] H.K.Jassal, J.S.Bagla, T. Padmanabhan, Mon.Not.Roy.Astron.Soc.Letters, 356, Lll (2005) 
|astro-ph/0404378| 



i q 



